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1*  Space  Charge  Limited  Flow  of  a  Thin  Electron  Beam  Confined  by 
a  Strong  Magnetic  Field 

An  approximate  analytic  theory  was  developed  and  implemented  numerically  for 
calculating  the  space  charge  limited  current  and  electric  field  of  a  thin  cylindrical 
beam  or  current  sheet  between  two  wide  parallel  electrodes.  The  flow  is  confined 
by  a  sufficiently  strong  magnetic  field.  Assuming  that  the  potential  and  current 
density  are  almost  homogeneous  in  the  direction  transversal  to  the  flow  we 
compute  the  beam  current  and  profile  by  a  variational  method.  We  find  that 
the  average  current  density  scales  as  the  reciprocal  of  the  beam  width  when  the 
latter  becomes  very  small.  The  total  cylindrical  beam  current  thus  decreases 
proportion! y  to  its  diameter  while  the  total  current  of  a  sheet  becomes  almost 
independent  of  the  width  in  this  regime. 

2*  Space  Charge  Limited  Flow  in  a  Rectangular  Geometry 

We  studied  the  spatial  structure  of  the  space  charge  limited  current  and  electric 
field  in  a  rectangle  of  arbitrary  aspect  ratio.  The  cathode  and  anode  form  two 
horizontal  sides  of  the  rectangle  and  a  strong  magnetic  field  forces  the  current  to 
flow  perpendicular  to  the  electrodes.  Using  conformal  mapping  techniques  we 
calculate  the  electric  field  outside  this  rectangle  for  any  given  potential  distribu¬ 
tion  on  its  vertical  boundaries.  Inside  the  current  rectangle  we  have  a  nonlinear 
Poisson  equation  with  extra  boundary  conditions  for  two  unknown!  functions: 
The  potential  and  the  current  density.  Both  exhibit  singular  behavior  at  the 
edges  of  the  rectangle,  A  semi  analytic  approximate  method  is  developed  for  this 
unusual  boundary  value  problem:  We  first  match  the  boundary  fields  inside  mid 
outside  the  current  region  and  then,  using  trial  functions  consistent  with  these 
matching  conditions,  wre  apply  the  least  square  technique  and  iterations  to  con¬ 
struct  the  solution  in  the  current  region.  The  analysis  of  the  flow  shows  that 
the  current  wings  are  similar  for  all  currents  wider  than  one  of  the  square  cross 
noindent.  There  is  also  evidence  that  the  total  current  does  not  vanish  when  the 
width  goes  to  zero.  The  method  of  calculation  appears  generalizable  to  various 
geometries  of  vacuum  and  solid  state  devices. 
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3.  Ionization  in  a  1-Dimensional  Dipole  Model 

We  study  the  evolution  of  a  one  dimensional  model  atom  with  ^-function  binding 
potential,  subjected  to  a  dipole  radiation  held  E(t)x  with  E(£)  a  2^/aj-periodic 
real- valued  function*  Starting  with  ip(x,  £  =  0)  an  initially  localized  state  and 
E(t)  a  trigonometric  polynomial,  complete  ionization  occurs;  the  probability 
of  finding  the  electron  in  any  fixed  region  goes  to  zero*  For  0)  compactly 
supported  and  general  periodic  fields,  we  construct  a  resonance  expansion.  Each 
resonance  is  given  explicitly  as  a  Gamow  vector,  and  2n/u  periodic  in  time 
and  behaves  like  the  exponentially  growing  Green’s  function  near  x  —  ±oo* 
The  remainder  is  given  by  an  asymptotic  power  series  in  t~ with  coefficients 
varying  with  x * 

4,  Excitation  and  Ionization  of  a  Simple  Two  Level  System  by  a 
Harmonic  Force 

A  simple  one-dimensional  quantum  system  with  two  attractive  ^-function  po¬ 
tentials  of  strength  q  at  x  ±  1  is  subjected  to  a  spatially  asymmetric  (as  in  the 
dipole  interaction)  harmonic  forcing  with  frequency  cj*  The  time  evolution  of 
the  system,  which  has  two  discrete  energy  levels  —w*  <  u)a  (depending  on  q) 
and  a  continuum  spectrum,  exhibits  a  rich  dynamics  including  regimes  where 
the  rate  of  ionization  becomes  very  small  due  to  the  ’inverse*  Ramsauer  effect 
in  electron-atom  collisions.  The  two-photon  ionization  with  u;  sb  ujs/2  can  be 
enhanced  when  uta  —  though  the  one-photon  ionization  is  not  affected  sig¬ 
nificantly  by  the  location  of  excited  level*  When  uj  is  very  close  to  the  oneor 
two-photon  resonance  the  ionization  rate  can  differ  greatly  from  that  given  by 
the  low  order  perturbation  theory  even  for  small  forcing  amplitude*  This  is 
caused  in  part  by  the  fact  that  the  dynamic  Stark  effect  has  a  strong  depen¬ 
dence  on  and  may  shift  the  resonance  frequencies  ujsi  uJfl  up  and  down.  When 
the  ground  state  decays  faster  than  the  excited  state  and  w  is  not  close  to 
the  excited  level  at  late  times  becomes  and  remains  more  populated  than  the 
ground  state.  The  occupations  of  the  bound  states  oscillate  with  a  frequency 
that  can  be  quite  lowcompared  with,  in  particular  in  the  case  w  w  uja  sb  w^/2, 
but  it  approaches  w  when  oj  when  uj  Our  analysis  is  based  on  the  an¬ 

alytic  structure  of  the  wavefunct ion’s  Laplace  transform  in  time.  It  considers 
one-  and  two- 7 photon’  processes;  the  higher  order  multiphoton  processes  can 
also  be  treated  by  our  computational  scheme  which  goes  beyond  the  low  order 
perturbation  theory. 

5*  Regularity  of  the  Laplace  Transform  and  Asymptotics  of  the  Wave 
Function  in  Time-Periodic  Potentials 

We  study  the  transition  to  the  continuum  of  an  initially  bound  quantum  par¬ 
ticle  in  Id,  d  —  1,2,3,  subjected,  for  £  >  0,  to  a  time  periodic  forcing  of  arbitrary 
magnitude.  The  analysis  is  carried  out  for  compactly  supported  potentials,  sat¬ 
isfying  certain  auxiliary  conditions.  It  provides  complete  analytic  information 
on  the  time  Laplace  transform  of  the  wave  function.  From  this,  comprehensive 
time  asymptotic  properties  (Borel  summable  transsenes)  follow* We  obtain  in 
particular  a  criterion  for  whether  the  wave  function  gets  fully  delocalized  (com¬ 
plete  ionization).  This  criterion  shows  that  complete  ionization  is  generic  and 
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provides  a  convenient  test  for  particular  cases.  When  satisfied  it  implies  absence 
of  discrete  spectrum  and  resonances  of  the  associated  Floquet  operator.  As  an 
illustration  we  show  that  the  parametric  harmonic  perturbation  of  a  potential 
chosen  to  be  any  nonzero  multiple  of  the  characteristic  function  of  a  measurable 
compact  set  has  this  property. 

0.  Transition  to  the  Continuum  of  a  Particle  in  Time- Periodic  Po¬ 
tentials 

We  present  new  results  for  the  transition  to  the  continuum  of  an  initially 
bound  quantum  particle  subject  to  a  harmonic  forcing.  Using  rigorous  expo¬ 
nential  asymptotics  methods  we  obtain  explicit  expressions,  as  generalized  Borel 
summable  transseries,  for  the  probability  of  localization  in  a  specified  spatial 
region  at  time  t.  The  transition  to  the  continuum  occurs  for  general  compactly 
supported  potentials  in  one  dimension  and  our  results  extend  easily  to  higher  di¬ 
mensional  systems  with  spherical  symmetry.  This  of  course  implies  the  absence 
of  discrete  spectrum  of  the  corresponding  Floquet  operator. 

7*  Propagation  Effects  on  the  Breakdown  of  a  Linear  Amplifier 
Model:  Complex- Mass  Schrodinger  Equation  Driven  by  the  Square 
of  Gaussian  Field 

Solutions  to  the  equation  dt£(x,t)  —  ^A£(tf,i)  —  AjS(;r, t)|2£{x, 0  are  inves¬ 
tigated,  where  £(x,  t)  is  a  complex  Gaussian  field  with  zero  mean  and  specified 
covariance,  and  m  /  (1  is  a  complex  mass  with  Im(m)  >  0.  For  real  m  this 
equation  describes  the  backscatterlng  of  a  smoothed  laser  beam  by  an  optically 
active  medium.  Assuming  that  S(x,  t)  is  the  sum  of  a  finite  number  of  indepen¬ 
dent  complex  Gaussian  random  variables,  we  obtain  an  expression  for  the  value 
of  A  at  which  the  g-tb  moment  of  |£(x,  £)|  wr.r.t.  the  Gaussian  field  S  diverges. 

This  value  is  found  to  be  less  or  equal  for  all  m  ^  0,  Im(m)  >  0  and  |m|  <  +oo 
than  for  |m|  =  H-oo,  i.e.  when  the  A£  term  is  absent.  Our  solution  is  based  on  a 
distributional  formulation  of  the  Feynman  path- integral  and  the  Paley- Wiener 
theorem. 

8.  Note  on  a  Diffraction- Amplification  Problem 

We  investigated  the  solution  of  the  equation  5* E(x^  t)—i DdlE[xt  £)  -  A|£(x,  t)|2E(x,  t ), 
for  x  in  a  circle  and  5{x,  t)  a  Gaussian  stochastic  field  with  a  covariance  of  a 
particular  form.  It  is  shown  that  the  coupling Ac  at  wdiieh  <  \E\  >  diverges  for 
t  >—  I  (in  suitable  units),  is  always  less  or  equal  for  D  >  0  than  D  =  0,  E. 

A.  Carlen,  M.  C.  Carvalho,  R,  Esposito,  J.  L.  Lebowitz  and  R.  Marra,  Droplet 
Minimizers  for  the  Cahn-Hilliard  Free  Energy  Function  We  proved  a  theorem 

characterizing  the  minimizers  in  a  model  for  condensation  based  on  the  Cahn 
Hilliard  free  energy  functional.  In  particular,  wre  exactly  determine  the  critical 
density  for  droplet  formation. 

9.  Phase  Transitions  in  Equilibrium  Systems:  Microscopic  Models 
and  Mesoscopic  Free  Energies 

We  describe  the  derivation  of  mesoscopic  free  energy  functionals  for  systems 
with  long-range  Kac  potentials  and  analyse  them  in  the  presence  of  phase  tran¬ 
sitions.  This  yields  information  about  the  arrangement  of  phases  in  some  binary 
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mixtures  as  well  as  explicit  criteria  for  the  stability  of  liquid  droplets  in  a  one- 
component  fluid  described  by  a  Cahn  Hilliard  type  of  free  energy  functional. 

10.  Ishig  Models  with  Long-Range  Antiferromagnetic  and  Short 
Range  Ferromagnetic  Interactions 

We  studied  the  ground  state  of  a  d  dimensional  Ising  model  with  both  long 
range  (dipole^like)  and  nearest  neighbor  ferromagnetic  (FM)  interactions.  The 
long  range  interaction  is  equal  to  r“p,  p  >  d,  while  the  FM  interaction  has 
strength  7,  If  p  >  d  +  1  and  J  is  large  enough  the  ground  state  is  FM,  while  if 
d  <  p  <  d+  1  the  FM  state  is  not  the  ground  state  for  any  choice  of  J.  In  d  =  1 
we  show  that  for  any  p  >  1  the  ground  state  has  a  series  of  transitions  from 
an  antiferromagnetic  state  of  period  2  to  2h-  periodic  states  of  blocks  of  sizes  h 
with  alternating  sign,  the  size  h  growing  when  the  FM  interaction  strength  J  is 
increased  (a  generalization  of  this  result  to  the  case  Q  <  p  <  1  is  also  discussed). 
In  d>  2  we  prove,  for  d  <  p  <  d+  1,  that  the  dominant  asymptotic  behavior  of 
the  ground  state  energy  agrees  for  large  J  with  that  obtained  from  a  periodic 
striped  state  conjectured  to  be  the  true  ground  state.  The  geometry  of  contours 
in  the  ground  state  is  discussed. 

11.  Sharp  Interface  Motion  of  a  Fluid  Binary  Mixture 

We  derive  hydrodynamic  equations  describing  the  evolution  of  a  binary  fluid 
segregated  into  two  regions,  each  rich  in  one  species, which  are  separated  (on 
the  macroscopic  scale)  by  a  sharp  interface.  Our  starting  point  is  a  Vlasov* 
Boltzmann  (VB)  equation  describing  the  evolution  of  the  one  particle  position 
and  velocity  distributions,  £),  %  —  1, 2.  The  solution  of  the  VB  equation  is 

developed  in  a  Hilbert  expansion  appropriate  for  this  system.  This  yields  incom¬ 
pressible  Navier-Stokes  equations  for  the  velocity  held  u  and  a  jump  boundary 
condition  for  the  pressure  across  the  interface.  The  interface,  in  turn,  moves 
with  a  velocity  given  by  the  normal  component  of  u . 

12.  Large  Deviations  for  Stochastic  Model  of  Heat  Flow 

We  investigate  a  one  dimensional  chain  of  2N  harmonic  oscillators  in  which 
neighboring  sites  have  their  energies  redistributed  randomly.  The  sites  —  N  and 
Ar  are  in  contact  with  thermal  reservoirs  at  different  temperature  and 
Kipnis,  Marehioro,  and  Presutti  [?]  proved  that  this  model  satisfies  Fourier’s 
law  and  that  in  the  hydrodynamical  scaling  limit,  when  N  -4  oo,  the  stationary 
state  has  a  linear  energy  density  profile  0(u),  u  €  [—1,1].  We  derive  the  large 
deviation  function  5(6>{u))  for  the  probability  of  finding,  in  the  stationary  state, 
a  profile  0(u)  different  from  &(u ).  The  function  S(9)  has  striking  similarities  to, 
but  also  large  differences  from,  the  corresponding  one  of  the  symmetric  exclusion 
process.  Like  the  latter  it  is  nonlocal  and  satisfies  a  variational  equation.  Unlike 
the  latter  it  is  not  convex  and  the  Gaussian  normal  fluctuations  are  enhanced 
rather  than  suppressed  compared  to  the  local  equilibrium  state.  We  also  briefly 
discuss  more  general  model  and  find  the  features  common  in  these  two  and  other 
models  whose  S(9)  is  known. 

13.  The  Asymmetric  Exclusion  Process  and  Brownian  Excursions 
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We  consider  the  totally  asymmetric  exclusion  process  (TASEP)  in  one  dimension 
in  its  maximal  current  phase.  We  show,  by  an  exact  calculation,  that  the  non- 
Gaussian  part  of  the  fluctuations  of  density  can  be  described  in  terms  of  the 
statistical  properties  of  a  Brownian  excursion.  Numerical  simulations  indicate 
that  the  description  in  terms  of  a  Brownian  excursion  remains  valid  for  more 
general  one  dimensional  driven  systems  in  their  maximal  current  phase* 

14.  Entropy  of  Open  Lattice  Systems 

We  investigate  the  behavior  of  the  Gibbs-Shannon  entropy  of  the  stationary 
nonequilibrium  measure  describing  a  one-dimensional  lattice  gas,  of  L  sites, 
with  symmetric  exclusion  dynamics  and  in  contact  with  particle  reservoirs  at 
different  densities.  In  the  hydrodynamic  scaling  limit,  L  to  infinity,  the  leading 
order  (O(L))  behavior  of  this  entropy  has  been  shown  by  Bahadoran  to  be  that 
of  a  product  measure  corresponding  to  strict  local  equilibrium;  we  compute  the 
first  correction,  which  is  0(1)*  The  computation  uses  a  formal  expansion  of 
the  entropy  in  terms  of  truncated  correlation  functions;  for  this  system  the  k-th 
such  correlation  is  shown  to  be  0(£_fc+1).  This  entropy  correction  depends 
only  on  the  scaled  truncated  pair  correlation,  which  describes  the  covariance 
of  the  density  field.  It  coincides,  in  the  large  L  limit,  with  the  corresponding 
correction  obtained  from  a  Gaussian  measure  with  the  same  covariance. 

15.  Fourier*s  Law  for  a  Harmonic  Crystal  with  Self-Consistent  Stochas¬ 
tic  Reservoirs 

We  consider  a  d-dimensional  harmonic  crystal  in  contact  with  a  stochastic 
Langevln  type  heat  bath  at  each  site.  The  temperatures  of  the  “exterior”  left 
and  right  heat  baths  are  at  specified  values  Tl  and  Tr,  respectively,  while 
the  temperatures  of  the  “interior”  baths  are  chosen  self- consistently  so  that 
there  is  no  average  flux  of  energy  between  them  and  the  system  in  the  steady 
state.  We  prove  that  this  requirement  uniquely  fixes  the  temperatures  and 
the  self  consistent  system  has  a  unique  steady  state.  For  the  infinite  system 
this  state  is  one  of  local  thermal  equilibrium.  The  corresponding  heat  current 
satisfies  Fourier’s  law  with  a  finite  positive  thermal  conductivity  which  can  also 
be  computed  using  the  Green-Kubo  formula.  For  the  harmonic  chain  (d  =  1) 
the  conductivity  agrees  with  the  expression  obtained  by  Bolsterli,  Rich  and 
Vissdier  in  1970  who  first  studied  this  model*  In  the  other  limit,  d  >  1,  the 
stationary  infinite  volume  heat  conductivity  behaves  as  (Zrfd)^1  where  A*  is  the 
coupling  to  the  intermediate  reservoirs*  We  also  analyze  the  effect  of  having  a 
non-uniform  distribution  of  the  heat  bath  couplings*  These  results  are  proven 
rigorously  by  controlling  the  behavior  of  the  correlations  in  the  thermodynamic 
limit* 

16.  Absolute  Continuity  of  Projected  SRB  Measures  of  Coupled 
Arnold  Cat  Map  Lattices 

We  study  a  d-dimensional  coupled  map  lattice  consisting  of  hyperbolic  toral  au¬ 
tomorphisms  (Arnold  cat  maps)  that  are  weakly  coupled  by  an  analytic  coupling 
map.  We  construct  the  Sinai- Ruelle-Bo  wen  measure  for  this  system  and  study 
its  marginals  on  the  tori.  We  prove  they  are  absolutely  continuous  with  respect 
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to  the  Lebesgue  measure  if  and  only  if  the  coupling  satisfies  a  nondegeneracy 
condition. 

17.  On  the  (Boltzmann)  Entropy  of  Nonequilibrium  Systems 

Boltzmann  defined  the  entropy  of  a  macroscopic  system  in  a  macrostate  M  as 
the  log  of  the  volume  of  phase  space  (number  of  microstates)  corresponding  to 
M.  This  agrees  with  the  thermodynamic  entropy  of  Clausius  when  M  spec¬ 
ifies  the  locally  conserved  quantities  of  a  system  in  local  thermal  equilibrium 
(LTE).  Here  we  discuss  Boltzmann’s  entropy,  involving  an  appropriate  choice  of 
macro- variables,  for  systems  not  in  LTE.  We  generalize  the  formulas  of  Boltz¬ 
mann  for  dilute  gases  and  of  Resibois  for  hard  sphere  fluids  and  show  that 
for  macro-variables  satisfying  any  deterministic  autonomous  evolution  equation 
arising  from  the  microscopic  dynamics  the  corresponding  Boltzmann  entropy 
must  satisfy  an  K- theorem. 

18,  A  Random  Matrix  Model  of  Relaxation 

We  consider  a  two  level  system,  <S2,  coupled  to  a  general  n  level  system,  Sn,  via 
a  random  matrix.  We  derive  an  integral  representation  for  the  mean  reduced 
density  matrix  p(t)  of  $2  in  the  limit  n  -¥  00,  and  we  identify  a  model  of 
$n  which  possesses  some  of  the  properties  expected  for  macroscopic  thermal 
reservoirs.  In  particular,  it  yields  the  Gibbs  form  for  p(oo).  We  consider  also  an 
analog  of  the  van  Hove  limit  and  obtain  a  master  equation  (Markov  dynamics) 
for  the  evolution  of  p(t)  on  an  appropriate  time  scale. 

ID,  Microscopic  Origins  of  Irreversible  Macroscopic  Behavior:  An 
overview 

Time-asymmetric  behavior  a s  embodied  in  the  second  law  of  thermodynamics 
is  observed  in  individual  macroscopic  systems .  It  can  be  understood  as  aris¬ 
ing  naturally  from  time-symmetric  microscopic  laws  when  account  is  taken  of 
a)  the  great  disparity  between  microscopic  and  macroscopic  scales,  b)  initial 
conditions,  and  c)  the  fact  that  what  we  observe  is  u typical”  behavior  of  real 
systems  not  all  imaginable  ones.  This  is  in  accord  with  the  ideas  of  Maxwell, 
Thomson,  and  Boltzmann  and  their  natural  quantum  extensions.  Common  al¬ 
ternate  explanations,  such  as  those  based  on  equating  irreversible  macroscopic 
behavior  with  the  ergo  die  or  mixing  properties  of  probability  distributions  (en¬ 
sembles)  already  present  for  chaotic  dynamical  systems  having  only  a  few  de¬ 
grees  of  freedom  or  on  the  impossibility  of  having  a  truly  isolated  system,  are 
either  unnecessary,  misguided  or  misleading.  Specie  features  of  macroscopic 
evolution,  such  as  the  diusian  equation,  do  however  depend  on  the  dynamical 
instability  (deterministic  chaos)  of  trajectories  of  isolated  macroscopic  systems. 
Time-asymmetric  behavior  as  embodied  in  the  second  law  of  thermodynamics 
is  observed  in  individual  macroscopic  systems.  It  can  be  understood  as  aris¬ 
ing  naturally  from  time-symmetric  microscopic  laws  when  account  is  taken  of 
a)  the  great  disparity  between  microscopic  and  macroscopic  scales,  b)  initial 
conditions,  and  c)  the  fact  that  what  we  observe  is  “typical”  behavior  of  real 
systems-  not  all  imaginable  ones.  This  is  in  accord  with  the  ideas  of  Maxwell, 
Thomson,  and  Boltzmann  and  their  natural  quantum  extensions.  Common  al¬ 
ternate  explanations,  such  as  those  based  on  equating  irreversible  macroscopic 
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behavior  with  the  ergodic  or  mixing  properties  of  probability  distributions  (en¬ 
sembles)  already  present  for  chaotic  dynamical  systems  having  only  a  few  de¬ 
grees  of  freedom  or  on  the  impossibility  of  having  a  truly  isolated  system,  are 
either  unnecessary,  misguided  or  misleading.  Specie  features  of  macroscopic 
evolution,  such  as  the  diusion  equation,  do  however  depend  on  the  dynamical 
instability  (deterministic  chaos)  of  trajectories  of  isolated  macroscopic  systems. 
20.  Product  Measure  Steady  States  of  Generalized  Zero  Range  Pro¬ 
cesses 

We  establish  necessary  and  sufficient  conditions  for  the  existence  of  factorizable 
steady  states  of  the  generalized  zero  range  process  on  a  periodic  or  infinite 
lattice*  This  process  allows  transitions  from  a  site  i  to  a  site  (i  +q)  involving 
(a  bounded  number  of)  multiple  particles  with  rates  depending  on  the  content 
of  the  site  i,  the  direction  q  of  movement,  and  the  number  of  particles  moving. 
We  also  show  the  sufficiency  of  a  similar  condition  for  the  continuous  time 
mass  transport  process,  where  the  mass  at  each  site  and  the  amount  transferred 
in  each  transition  are  continuous  variables;  we  conjecture  that  this  is  also  a 
necessary  condition. 

21*  Population  Dynamics  in  Spatially  Heterogeneous  Systems  with 
Drift:  The  Generalized  Contact  Process 

We  investigated  the  time  evolution  and  stationary  states  of  a  stochastic,  spa¬ 
tially  discrete,  population  model  (contact  process)  with  spatial  heterogeneity 
and  imposed  drift  (wind)  in  one-  and  two-dimensions.  We  consider  in  particu¬ 
lar  a  situation  in  which  space  is  divided  into  two  regions:  an  oasis  and  a  desert 
(low  and  high  death  rates).  Carrying  out  computer  simulations  we  find  that  the 
population  in  the  (quasi)  stationary  state  will  be  zero,  localized,  or  delocalized, 
depending  on  the  values  of  the  drift  and  other  parameters.  The  phase  diagram 
is  similar  to  that  obtained  by  Nelson  and  coworkers  from  a  deterministic,  spa¬ 
tially  continuous  model  of  a  bacterial  population  undergoing  convection  in  a 
heterogeneous  medium* 

22*  Behavior  of  Susceptible  Infected  Susceptible  Epidemics  on  Het¬ 
erogeneous  Networks  with  Saturation 

We  investigated  saturation  effects  in  susceptible-infected-susceptible  (STS)  mod¬ 
els  of  the  spread  of  epidemics  in  heterogeneous  populations*  The  structure  of 
interactions  in  the  population  is  represented  by  networks  with  connectivity  dis¬ 
tribution  P( ft), including  scale- free (SF)  networks  with  power  law  distributions 
P(k)  ~  fc’T  Considering  cases  where  the  transmission  of  infection  between 
nodes  depends  on  their  connectivity,  we  introduce  a  saturation  function  C(fe) 
which  reduces  the  infection  transmission  rate  A  across  an  edge  going  from  a 
node  with  high  connectivity  te.  A  mean  field  approximation  with  the  neglect 
of  degree-degree  correlation  then  leads  to  a  finite  threshold  \c  >  0  for  SF  net¬ 
works  with  2  <  7  <  3.  We  also  find,  in  this  approximation,  the  fraction  of 
infected  individuals  among  those  with  degree  k  for  A  close  to  A^.  We  investigate 
via  computer  simulation  the  contact  process  on  a  heterogeneous  regular  lattice 
and  compare  the  results  with  those  obtained  from  mean  field  theory  with  and 
without  neglect  of  degree-degree  correlations* 
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23-  Pair  Approximation  of  the  Stochastic  Susceptible- In  fee  ted- Recovered- 
Susceptible  Epidemic  Model  On  the  Hypercubic  Lattice 

We  investigated  the  time-evolution  and  steady  states  of  the  stochastic:  susceptible- 
infee  ted- recover  ed-susceptible(SIRS)  epidemic  model  on  one-  and  two-  dimen¬ 
sional  lattices.  We  compare  the  behavior  of  this  system,  obtained  from  computer 
simulations,  with  those  obtained  from  the  mean-field  approximation  (MFA)  and 
pair-approximation  (PA  )♦  The  former(latter)  approximates  higher  order  moments 
in  terms  of  first  (second)  order  ones.  We  find  that  the  PA  gives  consistently  bet¬ 
ter  results  than  the  MFA*  In  one  dimension  the  improvement  is  even  qualitative* 

24,  Large  Deviations  for  a  Point  Process  of  Bounded  Variability 

We  consider  a  one-dimensional  translation  invariant  point  process  of  density  one 
with  uniformly  bounded  variance  of  the  number  Nj  of  particles  in  any  interval 
/*  Despite  this  suppression  of  fluctuations  we  obtain  a  large  deviation  principle 
with  rate  function  F(p)  ~  —L~l  log  Prob(p)  for  observing  a  macroscopic  density 
profile  p(x),  xe[0,  Ij,  corresponding  to  the  coarse-grained  and  rescaled  density  of 
the  points  of  the  original  process  in  an  interval  of  length  L  in  the  limit  L  ->  oo. 

F(p)  is  not  convex  and  is  discontinuous  at  p  s  1,  the  typical  profile- 

25,  Point  Processes  with  Specified  Low  Order  Correlation 

Given  functions  p*(ri,  - .  r^),  defined  on  (Rd)k  or  (Zd)k  for  k  —  I,...,/!,  we 
wish  to  determine  whether  they  are  the  first  n  correlation  functions  of  some 
point  process  in  Rd  or  Zd ,  respectively,  and  if  so  what  we  can  say  about  the 
process.  We  give  partial  answers  to  these  questions  and  discuss  some  examples. 

26*  On  the  Realizability  of  Point  Processes  with  Specified  One  and 
Two  Particle  Densities 

We  investigated  and  gave  partial  answers  to  the  following  question:  given  one 
particle  and  pair  density,  pi(r%)  and  r^r^lR^  (or  does  there 

exist  a  point  process,  i*e.  a  probability  measure  on  points  in  (£d),  having 
these  densities? 

27.  On  the  Construction  of  Particle  Distribution  with  Specified  Single 
and  Pair  Densities 

We  discuss  necessary  conditions  for  the  existence  of  probability  distribution  on 
particle  configurations  in  d-dimensions  i*e.  a  point  process,  compatible  with 
a  specified  density  p  and  radial  distribution  function  g(r).  In  d  —  1  we  give 
necessary  and  sufficient  criteria  on  ptf(r)  for  the  existence  of  such  a  point  process 
of  renewal  (Markov)  type*  We  prove  that  these  conditions  are  satisfied  for  the 
case  <?(r)  -  0,  r  <  D  and  g(r)  —  l,r  >  D,  if  and  only  if  pD  <  e~l  :  the  maximum 
density  obtainable  from  diluting  a  Poisson  process.  We  then  describe  briefly 
necessary  and  sufficient  conditions,  valid  in  every  dimension,  for  pg(r)  to  specify 
a  determinantal  point  process  for  which  all  n-particle  densities,  pn(ri , rn), 
are  given  explicitly  as  determinants.  We  give  several  examples. 

28-  On  the  Distribution  of  the  Wave  Function  for  Systems  in  Thermal 
Equilibrium 
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For  a  quantum  system,  a  density  matrix  p  that  is  not  pure  can  arise,  via  averag¬ 
ing,  from  a  distribution  p  of  its  wave  function,  a  normalized  vector  belonging  to 
its  Hilbert  space  Jf.  While  p  itself  does  not  determine  a  unique  p,  additional 
facts,  such  as  that  the  system  has  come  to  thermal  equilibrium,  might.  It  is  thus 
not  unreasonable  to  ask,  which  p,  if  any,  corresponds  to  a  given  thermodynamic 
ensemble?  To  answer  this  question  we  construct,  for  any  given  density  matrix 
p ,  a  natural  measure  on  the  unit  sphere  in  Jf,  denoted  GAPp ,  We  do  this 
using  a  suitable  projection  of  the  Gaussian  measure  on  Jf  with  covariance  p . 
We  establish  some  nice  properties  of  GAPp  and  show  that  this  measure  arises 
naturally  when  considering  macroscopic  systems.  In  particular,  we  argue  that  it 
is  the  most  appropriate  choice  for  systems  in  thermal  equilibrium,  described  by 
the  canonical  ensemble  density  matrix  =  {l/Z)  exp(— flH).  GAPp  may  also 
be  relevant  to  quantum  chaos  and  to  the  stochastic  evolution  of  open  quantum 
systems,  where  distributions  on  Jf  are  often  used. 

29.  Canonical  Typicality 

It  is  well  known  that  a  system  5  weakly  coupled  to  a  heat  bath  B  is  described 
by  the  canonical  ensemble  when  the  composite  S  B  is  described  by  the  mi- 
crocanonical  ensemble  corresponding  to  a  suitable  energy  shell.  This  is  true 
for  both  classical  distributions  on  the  phase  space  and  quantum  density  matri¬ 
ces.  Here  we  show  that  a  much  stronger  statement  holds  for  quantum  systems. 
Even  if  the  state  of  the  composite  corresponds  to  a  single  wave  function  rather 
than  a  mixture,  the  reduced  density  matrix  of  the  system  is  canonical,  for  the 
overwhelming  majority  of  wave  functions  in  the  subspace  corresponding  to  the 
energy  interval  encompassed  by  the  mic mean oni cal  ensemble.  This  clarifies, 
expands,  and  justifies  remarks  made  by  Schrodinger  in  1952. 

30.  Using  Kinetic  Monte  Carlo  Simulations  to  Study  Phase  Separa¬ 
tion  in  Alloys 

We  review  recent  extensions  of  the  kinetic  Ising  model  used  to  investigate  phase 
separation  in  binary  alloys.  Firstly,  vacancies  are  included  to  model  the  diffusion 
of  the  atoms  on  the  microscopic  scale  more  real  is  tic  ally.  These  can  change 
the  coarsening  rate  and  the  coarsening  mechanism.  Secondly,  the  lattice  is 
allowed  to  deform  owing  to  the  different  sizes  of  the  atoms  and  the  resulting 
misfit  between  precipitates  and  matrix.  The  deform  ability  of  the  lattice  induces 
long-range  elastic  interactions  between  the  atoms.  These  change  the  shape, 
orientation,  and  arrangement  of  the  precipitates.  The  growth  of  the  precipitates 
need  not  follow  the  H(£)a  law. 

31.  Stability  of  Solutions  of  Hydrodynamic  Equations  Describing  the 
Scaling  Limit  of  a  Massive  Piston  in  an  Ideal  Gas 

We  analyze  the  stability  of  stationary  solutions  of  a  singular  Vlasov  type  hydro- 
dynamic  equation  (HE).  This  equation  was  derived  (under  suitable  assumptions) 
as  the  hydrodynamical  scaling  limit  of  the  Hamiltonian  evolution  of  a  system 
consisting  of  a  massive  piston  immersed  in  an  ideal  gas  of  point  particles  in  a 
box.  We  find  explicit  criteria  for  global  stability  as  well  as  a  class  of  solutions 
which  are  linearly  unstable  for  a  dense  set  of  parameter  values.  We  present 
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evidence  (but  no  proof)  that  when  the  mechanical  system  has  initial  conditions 
‘'close”  to  stationary  stable  solutions  of  the  HE  then  it  stays  close  to  these  so¬ 
lutions  for  a  time  which  is  long  compared  to  that  for  which  the  equations  have 
been  derived.  On  the  other  hand  if  the  initial  state  of  the  particle  system  is 
close  to  an  unstable  stationary  solutions  of  the  HE  the  mechanical  motion  fol¬ 
lows  for  an  extended  time  a  perturbed  solution  of  that  equation:  we  find  such 
approximate  periodic  solutions  that  are  linearly  stable. 

32,  Percolation  in  the  Harmonic  Crystal  and  Voter  Model  in  Three 
Dimensions 

We  investigate  the  site  percolation  transition  in  two  strongly  correlated  systems 
in  three  dimensions:  the  massless  harmonic  crystal  and  the  voter  model.  In  the 
first  case  we  start  with  a  Gibbs  measure  for  the  potential,  U  =  -  ^<xy>  {$(*)  —  ^(y))2, 
x,  y  e  Z3,  J  >  0  and  0(x)  €  R,  a  scalar  height  variable,  and  define  occupation 
variables  ph(z)  =  1,(0)  for  $(%)  >  h{<  h ).  The  probability  p  of  a  site  being 
occupied,  is  then  a  function  of  h .  In  the  voter  model  we  consider  the  stationary 
measure,  in  which  each  site  is  either  occupied  or  empty,  with  probability  p .  In 
both  cases  the  truncated  pair  correlation  of  the  occupation  variables,  G(x  —  y), 
decays  asymptotically  like  \x  ~  y|_1.  Using  some  novel  Monte  Carlo  simulation 
methods  and  finite  size  scaling  we  find  accurate  values  of  pc  as  well  as  the  critical 
exponents  for  these  systems.  The  latter  are  different  from  that  of  independent 
percolation  in  d  —  3,  as  expected  from  the  work  of  Weinrib  and  Halperin  [WH] 
for  the  percolation  transition  of  systems  with  G(r)  ^  r~a  [A.  Weinrib  and  B. 
Halperin,  Phys.  Rev.  B  27,  413  (1983)].  In  particular  the  correlation  length 
exponent  v  is  very  close  to  the  predicted  value  of  2  supporting  the  conjecture 
bv  WH  that  v  —  -  is  exact. 

■  a 

33,  Percolation  Phenomena  in  Low  and  High  Density  Systems 

We  consider  the  2D  quenched-disordered  q-  state  Potts  ferromagnets  and  show 
that  at  seif-dual  points  any  amalgamation  of  q  —  I  species  will  fail  to  percolate 
despite  an  overall  (high)  density  of  1  —  g”1.  Further,  in  the  dilute  bond  version 
of  these  systems,  if  the  system  is  just  above  threshold,  then  throughout  the  low 
temperature  phase  there  is  percolation  of  a  single  species  despite  a  correspond¬ 
ingly  small  density.  Finally,  we  demonstrate  both  phenomena  in  a  single  model 
by  considering  a  “perturbation”  of  the  dilute  model  that  has  a  self  dual  point. 

We  also  demonstrate  that  these  phenomena  occur,  by  a  similar  mechanism,  in 
a  simple  coloring  model  invented  by  0.  Haggstrom. 

34,  Space  Charge  Limited  Two-dimensional  Unmagnet i^ed  Flow  in  a 
Wedge  Geometry 

This  paper  studies  the  space  charge  limited  current  in  an  infinite  wedge  geometry 
in  two  dimensions.  This  geometry  permits  a  reduction  of  the  problem  to  a 
set  of  easily  solved  ordinary  differential  equations-  The  system,  though  very 
simplified,  exhibits  features  similar  to  those  expected  to  occur  in  many  realistic 
systems  with  inhomogeneous  electric  fields.  We  obtain,  in  particular,  a  universal 
form  for  the  particle  trajectories  and  a  non-monotone  charge  distribution  with 
accumulation  at  both  the  cathode  and  the  anode.  The  explicit  solution  of  the 
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model  can  be  useful  for  testing  numerical  schemes.  The  case  of  a  very  low  density 
current  is  also  considered.  Relaxation  of  the  geometrical  limitations  of  the 
model  are  studied  using  conformal  mapping  techniques.  Possible  applications 
to  realistic  systems,  which  can  be  tested  by  simple  experiments,  are  presented. 
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